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Abstract 

Phase separation in binary and ternary fluids is studied using a two dimensional 
Lattice Gas Automata. The lengths, given by the the first zero crossing point of 
the correlation function and the total interface length is shown to exhibit power law 
dependence on time. In binary mixtures, our data clearly indicate the existence of 
a regime having more than one length scale where the coarsening process proceeds 
through the rupture and reassociation of domains. In ternary fluids; in the case of 
symmetric mixtures there exists a regime with a single length scale having dynamic 
exponent 1/2, while in asymmetric mixtures our data establish the break down of 
scale invariance. 
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I. INTRODUCTION 



A mixture of fluids phase separate into domains when quenched below its critical tem- 
perature. Since many systems of scientific and technological interest are multi component 
mixtures, the phase behavior of fluid mixtures are of current interest. The equilibrium state 
of incompatible fluids is one in which the pure phases are separated by a single connected 
interface. However, in the thermodynamic limit, starting from a mixed phase this equilib- 
rium is never reached. In view of this, the kinetics of the phase separation process gains 
importance. 

In the case of phase separating binary fluids, experiments (1; 2; 3; 4) and numerical 
simulations (5; 6; 7) have clearly established the importance of hydrodynamics in the de- 
termination of late time domain growth laws. In spite of all efforts, a complete theoretical 
understanding of this highly nonlinear phenomena remains unsolved (8). This is mainly due 
to the fact that studying kinetics of liquid phase separation involves the solution of mutually 
coupled equations; the Navier-Stokes equation for the flow and the equation of continuity 
for the order parameter. 

Early experimental and theoretical studies of these binary fluid systems, assume that 
there exists, at late times, a dynamical scaling regime exhibiting similar behavior under an 
appropriate rescaling of time and length scales. Dynamical scaling is characterized by the 
single time-dependent length scale R(t). The domain growth follows a simple and generic 
algebraic form, R(t) oc t a ,where a represents the exponent characteristic of the universality 
class to which the system belongs. Scaling and dimensional analysis by Siggia (9) , Furukawa 
(10), San Miguel (11) and more recently by Bray (8) addresses this question of the growth 
exponent taking this length scale R(t) to be the average size of the ordering domains. In two 
dimensions, in the case of minority B phase separating from a A rich mixture , the domains 
coarsen via evaporation-condensation process leading to the exponent a = 1/3. The late 
time growth is governed by droplet coalescence leading to the exponent a — 1/2. In the case 
of symmetric binary mixtures, thermal fluctuations can drive the initial a = 1/3 regime to 
a a = 1/2 regime (11). Most of the studies, however show a cross over from the a = 1/3 to 
an inertial a = 2/3 regime predicted by Furukawa (10). 

More recently many researchers have pointed out the absence of scaling in two dimensional 
phase separating fluids. It is shown that competition between diffusive and hydrodynamic 
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growth leads to breakdown of scale-invariance in symmetric binary fluids (15). In the 
viscous hydrodynamic regime, scaling is observed only in the case of coarsening through 
coalescence. Even after starting from a droplet state one can enter a bicontinuous state by 
coalescence induced coalescence mechanism (16) leading to a breakdown of scaling. Alter- 
natively starting from a droplet state scaling is maintained, in symmetric binary mixtures, 
if the droplet morphology is self sustaining (18). In the inertial hydrodynamics regime ,full 
scaling is recovered (17; 18). Excellent scaling is observed also in the cross over regime from 
the viscous hydrodynamic to the inertial hydrodynamic regime (19). 

Variety of techniques have been used to simulate growth kinetics in the binary immiscible 
fluids. Direct simulations of model H have been used to study the domain growth and 
scaling using variously defined length scales (8; 17). Lattice Boltzmann simulations have 
been particularly useful in exploring the late time hydrodynamic regime (12; 13; 14). One 
problem with this technique is that it does not include thermal fluctuations. On the other 
hand thermal fluctuations are inherent in Lattice Gas models. Rothman and Keller proposed 
a Lattice Gas Model (RK model) for the binary immiscible fluids (20). This model has been 
used for simulating 2-D binary fluid phase separation at different overall fluid densities 
(21; 22). 

Though most of the techniques mentioned in the previous paragraph have been extended 
to study a mixture of two fluids and a surfactant (22; 23; 24), our understanding of phase 
separation wherein all three are fluid components is rather limited. Molecular Dynamics 
simulations by Laradji et.al (25) seem to indicate that hydrodynamic flow is not likely 
to control the separation process in ternary fluid system. This study indicates that the 
ternary system at late times reaches a dynamical scaling regime during which the domains 
show a growth law R(t) oc t 1 ^ 3 , in agreement with the classical theory of Lifshitz and 
Slyozov. Gunstensen and Rothman extended the Lattice Gas Model (GR Model) to study 
ternary immiscible fluids (LGTF) (27) but do not make any comments about the dynamical 
exponents. In this model the total energy function is the sum of the work done by the "color 

— * — * — * — * 

flux" fi of each component against its "color field" Qi and is given by W = J^^ifi ' Qi- 
More recently, a level set method was proposed for the study of phase separation in fluid 
mixtures (26). In this method one assumes convective terms to dominate the coarsening; 
local volume fraction function <f>i(x, t) of the components is then coupled to the local velocity 
field v through the kinetic equation, 
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(1) 



The fluid velocity v satisfies Navier-Stokes equation 
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Here P is the pressure ,u is the kinematic viscosity and p is the density. The interfacial 
energy between the domains enters the equations only through the external force F; 



i=A,B,C 

where n(fa) is the curvature of the fa domain interface, a is the surface tension defined 
such that for the AB interface the surface tension is a\ + a l B . 

The zero contour of the function fa specifies the interface of the domains. Immiscible three 
component mixtures with majority components A and B having the same volume fraction 
and a minority C phase was studied using this method (26). They find a power law 
dependence of interface length with dynamic exponents in the range of 0.5 ~ 0.6. 

In this paper we discuss Lattice gas simulations of the binary and ternary mixtures phase 
separating in two dimensions. We calculate cluster size distribution, total interface length 
and the density density correlation functions for each of the components. We see that, in 
some cases, the interface length and the first zero of the correlation function; though shows 
power law dependence on time, do not have the same dynamic exponents, thus violating 
scaling. 

In binary mixtures large difference in the parameters Oi of the two components can lead to 
break down of scaling with the domains coarsening through rupture and coalescence. In the 
symmetric ternary mixtures with equal <7j and equal volume fractions, we see the average size 
of the ordering domains, R(t) oc t 1 ^ in agreement with the results reported earlier. At late 
times, we observe a crossover to a t 1 ^ 2 regime which is not reported in previous simulations. 
As we reduce the volume fraction of one of the components this exponent changes to 2/3, 
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consistent with the inertial regime in two component mixtures. In the asymmetric case; 
wherein one of the component is miscible in the other two, we observe a R(t) oc t 1 / 2 and a 
cross over to t 2 l z as the volume fraction of the solute is reduced to zero. 

By choosing a substrate which has minimum interaction with the fluid one should be 
able to conduct experiments on fluid separation problem in two dimensions. Such quasi two 
dimensional geometries have been used for studying three fluid phase separation (4). 

The paper is organized as follows: In section I, we introduce the problem and summarize 
previous results. A detailed account of the Lattice Gas Automaton model and the algorithm 
for computing the cluster size and interface length are given in section II. Next, we discuss 
the results of our simulations on the binary fluids. Section IV deals with the symmetric 
and asymmetric ternary fluids. In section V,we conclude the paper with a summary of the 
results. 

II. THE MODEL 

Lattice Gas Automaton (LGA) is an alternative numerical description of fluid flow dy- 
namics. This model approximates reality by constraining motions and collisions of fluid 
particles to a lattice, where each particle represents a finite mass of fluid. Lattice-gas au- 
tomaton models represent, in the continuum, the incompressible Navier-Stokes equation 
correctly and are extensively used in simulations of fluids. 

In this paper we discuss simulations using the FHP model for fluids. Much has been 
written about the FHP model which was first introduced in 1986 by Frisch, Hasslacher 
and Pomeau (28). The fully discrete microscopic dynamics of the FHP model maps into 
the macroscopic behavior of hydrodynamics. A modified version of the FHP model for 
immiscible fluid was introduced by Rothman and Keller (20). In this model, which is defined 
on a hexagonal lattice, particles at the lattice sites are allowed to take 7 possible velocities, 
cl,i G {0, 1,2,3,4,5,6} . Boolean variables indicate the occupation number at a particular 
site at a given time. The dynamics is such that no more than one particle enters the same 
site at the same time with the same velocity; the exclusion principle. The particles undergo 
collision that conserves, at every site, the total number of particles, their total momentum 
and the total kinetic energy after each time step. 

Multi phase flows in which different species of fluids coexist, move and interact are an 
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Fig. 1 The state of a lattice point is given by the vectors <fi = {1 5 0, 1, 0, 0, 1, 0} , 
q{ = {0, 0, 1, 0, 0, 0, 0}, q-i = {1, 0, 0, 0, 0, 0, 0} and q = {0, 0, 0, 0, 0, 1, 0} 



important domain of application of the lattice gas approach. Models for miscible, immiscible 
and reactive flows have been proposed. Interaction between the species are the key ingre- 
dient of multi-phase flows and depending on the nature of the physical processes, different 
interactions will be considered. 

The Lattice gas approach has the important advantage that the interface between the 
different fluids appear naturally as a consequence of the way the fluids are modeled in terms 
of particles. In this paper, we represent the different phases by assigning the particles a 
"color" . Now the state of any vertex at r is specified by a total density vector 0(f) and the 
color vectors qi{r) with Boolean variables 4>j an d q%j as the components of this vector. In 
the case of ternary mixtures, % can take 3 values +1,-1 and . This is illustrated clearly 
in figure 1. We define the average density as the ratio of total number of particles to the 
maximum number of particles possible. This is given by d — ^lZ*oj=i < f > i( r j)/('fN), where 
iV is the total number of sites in the lattice. 

Drawing analogies from electrodynamics, we define a "color flux" and "color field" at 
every lattice point for each type of particles. In the case of ternary mixture the color flux 
at a lattice site is given by 




(4) 



3=0 




(5) 



3=0 
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And the local color gradients or fields are defined to be 



j=5 k=6 



/i( r ) = Yl ^ Yl ^(r + Cj) 



(7) 



j=0 k=0 



j=5 k=6 



(8) 



j=0 k=0 



j=5 k=6 

fo( r ) = J2^Y1 9o*(r + Cj) (9) 

j=0 k=0 

Following the Gunstensen and Rothman (GR model), we write the work done by the flux 
against the field to be: 



where a, are the parameters which determine the surface tension between different phase 
boundaries. The three phases are labeled by % — ±1, 0. 

The simulation proceeds in two steps, Collision and Translation. During the collision 
process, particles at a given lattice site can exchange their velocity and color. The collision 
rules are such that the work performed by the flux against the field is minimized , subject to 
the constraints of the exclusion principle as well as the principles of mass , momentum and 
particle conservation. There can be a number of configurations which satisfy the above con- 
straints, and in principle one should pick up that configuration which minimizes energy from 
this full set of possible configurations. In practice only a "limited" number of configurations 
are considered. In our simulations we take into account the following possibilities. 

If a site has more than three particles or if their total momentum is non zero, then we 
pick up all possible pairs and an attempt is made to exchange their colors. The exchange is 
accepted only if it lowers the energy. 

When 2 particles occupy the same site with opposite velocities there are two cases to 
consider. In one both the particles are of the same color. Then a rotation of the configuration 
by 0, 60 or 120 degrees is attempted as shown in figure 2. The other possibility is that the 



W = -(aJ^.Qiir) + a-if-W.Ql^r) + a f (r).Q (r)) 



(10) 
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Fig. 2 Possible configurations in a two particle collision 

two are of different colors. Then we need to check the possibilities for 0, ±60, ±120 and 180 
degree rotation. 

When three particles collide at an angle of 120 degrees there are three possibilities. 
When all the particles have the same color, they could just reverse their directions or retain 
their velocities. When the colors are not the same, similar to two particle collisions, the 
configuration could now be rotated by 0, ±60, ±120 and 180 . In all the above cases the 
lowest energy configuration is chosen as the final state. 

There is one special case of three particle collision in which two particles have opposite 
velocities and the third one is a zero velocity particle. In this case, apart from rotating the 
configuration an exchange of color between pairs is also attempted. 

Next step in the simulation is the translation in which all particles are moved in the 
direction of their velocities. This collision and translation process completes one time step. 

To understand the various competing terms in this model we first look at a two component 
fluid mixture and introduce the variables qk = qik — Q-ik, <Pk = qik + Q-ik, in the direction 
k. The two components are now designated by qk — ±1. Substituting these variables in the 
equations for f\, /_i and Qi, Q-i, we can obtain from equation 10 the work done to be 

W = -(<7 1 / 1 Q 1 + (7_ 1 /_ 1 Q_ 1 ) 

= ~(y] c fc <?fc I (oi + a-i) X Ci ^ qi(r + + (<n - <7_i) ^ Q ^ r + Q ) 

k \ i l i I 

+ Ck( f )k [ ( ai + <7 - 1 ) Ci ^ r + C *) ~ ( <Tl ~ <7 - 1 ) S Ci qi ( r + Q ) ) X 11 ) 

k \ i I i I / 

Momentum conservation implies that the last two terms in the above equation do not 
contribute to the dynamics. The first term decides the energy of the interface between the 
two phases. The second term favors movement of particles toward higher or lower densities 



depending on their color. For example when o\ — o~\ > this terms imply higher diffusivity 
of particles with q = — 1. Thus increasing \a\ — cr_i| is like increasing the "temperature". 

The earlier lattice gas model for two component fluids by Rothman and Keller (RK 
model) used <j\ = <7_i = a. In this model the work done is given by 

W = -of- Q, (12) 

where Q(r) = J^jZlcjQj and f( r ) = J2jZl c j J2k=o Qk(r + Cj). We can then intro- 
duce thermal fluctuation by defining an inverse temperature like parameter f3 such that 
the new configuration at every time step is accepted with a probability proportional to 
exp — (3 {Woid — Wnew) (22). Unless specified otherwise the results discussed in this paper 
are obtained using the GR model. 



A. Correlation functions and domain distribution functions 

We define the pair correlation functions 

\ k k I \ k I \ k 

Cj(f,0 = (h^hjix + ^t)) - (h t {x)} {h 3 {x)) , (14) 

Here the subscripts i,j represent the color of the particles with i,j = ±1,0. Field values 
hi determines the phase at every lattice points. A lattice point belongs to the i th phase if hi 
has a value higher than the other two. This means that most of the particles there belong 
to the i th color. We calculate the length Rij(t) as the first zero of the correlation function 

r~fh 

For the immiscible case wherein all three phases have the same surface tension, to calcu- 
late the border length and cluster size of the domains of a particular phase, say A, we use 
the following algorithm. 

All lattice sites are picked sequentially. If a particular site belongs to the phase A and is 
not assigned a cluster number we give it a new number. All the neighbors of this site are 
then checked. If the neighbor belongs to the phase A, three possibilities arise. In one, the 
site does not have an assigned cluster number. We will then give this neighbor the same 
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Fig. 3 (top) Rules for determining the border. Square represent C phase, Open circles 
represent B phase and the hashed circles represent A phase. The dark squares and circles are the 
border sites between the A-C phase and the B-C phase, (bottom) Snap shot from the simulation 
of ternary fluid with one component miscible in the other two. The dark regions are the border 

between th A-C (white) and B-C (gray) domains. 

cluster number as the present site. In the second case, the neighbor has a number which 
is the same as the site . Here we do not have to do anything more. In the third case the 
neighbor has a number which is different from that of the present site. That is, the neighbor 
belongs to a different cluster. In that case, the clusters are merged by assigning smaller of 
the two numbers to all sites in the two clusters. A lattice site which belong to a particular 
phase is a border site if one of its neighbors is of a different phase. 

In the case, wherein one of the phases is equally miscible in the other two; for example let 
us consider C miscible in A and B, we have two phases, the A-C phase and the B-C phase, 
see figure 3. To label the domains we then follow the same procedure as before. Any site 
belongs to the A-C phase if it is an A site or a C site whose neighbors are all A or C. For 
example region c in figure 3 belongs to the A-C phase. The border of the A-C domain could 
be a A or C site. If it is an A site, it should have at least one B neighbor (region d in figure 
3) or a C neighbor with the next nearest neighbor as B (region a in figure 3). If the border 
site belongs to the C phase then it should either have a neighbor which is a B site (region 
b in figure 3) or a C neighbor with the next nearest neighbor as B (region a in figure 3). 
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Fig. 4 The first zero of the correlation function R (open circles) and the total border length C 
(filled squares) for a binary fluid, with g\ = cr_i = 1, as a function of time. The average density 
of the fluids is d = 0.55. The curves are shifted to make the comparison easy. Continuous lines 
with slope 2/3 and 1/3 are given as guide to the eye 

In the simulations, to prepare the initial configuration a particle is placed at a randomly 
chosen site. The velocity of the particle is chosen randomly from the seven possible directions 
such that the total momentum summed over all particles is zero. The color of the particle 
is again assigned randomly depending on the volume fraction of each component. Once 
the initial configuration is prepared the ensuing dynamics conserved the total number of 
particles, the net momentum, the total kinetic energy and the volume fraction of each 
components. 

Simulations are performed on a triangular lattice with 360000 lattice points with periodic 
bondary conditions. The results are averaged over 10 initial configurations. The dynamic 
exponents are obtained by a power law fit over one decade in time. 

III. TWO COMPONENT FLUIDS 

In this section we reexamine the well studied case of binary fluid coarsening following a 
critical quench. The fluids are labeled by i — ±1. Let us first look at the case U\ = <r_i = 1. 
The time dependence of average domain size R(t) which scales the same way as R^ is shown 
in figure 4. At early time we see Lifshitz-Slyozov-Wagner(LSW) t 1//3 (8) growth. At late 
times the growth of domains is controlled by the inertial hydrodynamics and we have the 
well known t 2 ^ growth law (8). 

To check for scaling we compare the length R given by the first zero of the correlation 
function with the total border length C(t) . As can be seen from figure 4, at late time, in 
the inertial hydrodynamic regime, the length C also exhibits power law dependence with 
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Fig. 5 Dynamic exponents of ctR and ac as a function of cr_i for a\ = 1. The open circles 
represent ac and the dark squares represent an 




Fig. 6 Time evolution of binary mixture for d=0.55, u\ = 1.0, cr_i = —0.4. The snap shots are 

taken at (a)* = 25118 and (b)t = 39810. 

a dynamic exponent 2/3. This agrees with earlier studies on binary fluid phase separation 
which predict scaling to hold in the inertial regime (17; 18). 

We now investigate the more general case of o\ ^ (J-\ with cr_! < <j\. The dynamic 
exponents cxr and ac characterising the time dependence of R and C seem to vary differently 
with a i — <J-\ This is depicted in figure 5. This results imply break down of scaling . The 
absence of scaling is clearly manifested in the snap shots given in figure 6, the coarsening 
process is through the breaking and joining of domains. Since the length R(t) is decided 
by the big clusters it has a dynamical exponent characteristic of coalescence. On the other 
hand the total border length C(t) has significant contribution from the small clusters which 
coarsen through evaporation condensation resulting in an exponent 1/3. 




Fig. 7 Time evolution of binary mixture in the RK model for d=0.55, a = 1.0, /3 = 0.03. The 
snap shots are taken at (a)t = 25118 and (b)t = 39810. 
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Fig. 8 Dynamic exponents of an and ac as a function of (3 for a = 1 in the RK model. The 
open circles represent ac and the dark squares represent 



Similar change in the dynamical exponents of R was observed as a function of the inverse 
temperature like parameter (3 in the RK model (22). For low (3 values the coarsening 
mechanism very similar to that operating in the case of U\ ^ o\_iin the GR model. This 
is shown in figure 7. The time dependence of R and C shown in figure 8 confirms absence 
of scaling. These exponents indicate the presence of more that one coarsening mechanism 
operating at the late time regime. 

In the case of two fluid mixtures with unequal volume fraction we see the initial LSW 
regime crossing over to the a — 1/2 regime controlled by the droplet coalescence mechanism. 
The first zero of the correlation function R and the total border length C in this regime seems 
to indicate existence of simple scaling (18). 



IV. THREE COMPONENT FLUIDS 



In this section, we look at two separate cases of phase separating three component fluids. In 
one, we consider a symmetric ternary mixture with all the components having equal volume 
fraction and with the parameters <7j in equation-10 set to o\ = <7_i = a . This amounts to 
having the same surface tension for the A — B, B — C and A — C interfaces. The other case 
of interest is the asymmetric mixture, where we choose the parameters to be o\ = U-\ and 
(j < 2(<7i + <r_i). This choice of parameters ensures that one of the components (solute) is 
equally soluble in the other two. We will now discuss these two cases in detail. 
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Fig. 10 Correlation length as a function of time in the immiscible three component case for 
densities are 0.45,0.5 and 0.55. a = .5 line is given as a guide to the eye. 



1. Symmetric case 

Starting from the mixed state, the symmetric ternary mixture, when quenched below the 
transition temperature, evolves to form droplets of individual components. Snap shots from 
the simulation, shown in figure 9, clearly establish the existance of sharp interfaces between 
the components. The first zero of correlation function R and the interface length C are 
plotted in figure 10 as a function of time. 

Initial regime corresponds to the formation of droplets. In figure-10 we see a region with 
R(t) ~ t 1//2 , where this droplets coalesce to form bigger domains. If the late time growth is 
due to the coalescence by droplet diffusion, as it is believed to be, then the following simple 
argument shows that the exponent should be 1/2. 

If R is the typical radius of a droplet, we have the droplet number density n ~ v/R 2 , 
where v is the volume fraction of the components in two dimensions. The time for a droplet 
to diffuse a distance of order of its radius is t R = R 2 /D , where D is the diffusion coefficient. 
The area swept out by the drop in time t (for t > t R )is of the order R 2 t/t R ~ Dt . If t c 
is the coalescence time , the expected number of drops in an area Dt c is of the order unity, 
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Fig. 11 The first zero of the correlation function R(open circles) and the total border length C 
(filled squares) for a symmetric ternary fluid with o\ = <r_i = &q as a function of time. Average 
density of the fluids is d=0.55. Continuous lines with slopes 1/3 and 1/2 are given as guide to the 

eye 



which means nDt c = 1. Since in two dimensions, the diffusion coefficient, D ~ ksT /rj, does 
not have any dependence on the droplet size, 



R 2 r)R 

vD vk B T K J 

This implies that R grows with time as R ~ (i^rt) 1 ^ (9). 

Another mechanism for droplet coalescence is interfacial diffusion. At non zero temper- 
atures the domains fluctuate from their circular shape. The area explored by the domain 
in time t goes as D b t, where D b is the diffusion coefficient for the interface. The coalescence 
time t c is then given by 



here the length R is such that tcR^ = area of the droplet. Since the interface diffusion 
coefficient is not a function of R , we get L(t) ~ t 1 / 2 . Thus both the droplet diffusion and 
interfacial diffusion could give the same dynamical exponent. 

To explore the possibility of the second mechanism we studied the stability of a drop to 
fluctuations from the circular shape. As shown in the appendix, in two dimensions, surface 
fluctuations die out exponentially, both in the case of a drop and a strip (11) of one fluid in 
another. This implies that in the a — 1/2 regime, the dominant mechanism for coarsening 
is the droplet coalescence by diffusion. 
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Fig. 12 Time evolution of three fluid mixture for d=0.55 . (a) and (b) are the snap shots from 
the asymmetric ternary mixture phase separation at t = 25118 and t = 39810 with <j\ = a_i = 1 
and do = —1.5. Note the breaking and reorganization of the domains with time. A-C region is 
shown in black and B-C region in white. For comparison we show snap shots (c) and (d) from the 
phase separation of a binary mixture having the same density at the same time steps with 

(71 = = 1. 

Droplet coalescence is the only mechanism operating here is further confirmed by the 
existence of scaling. In figure 11 we compare the time dependence of the total border length 
Ca of the domains of component A with that of the first zero of the correlation function 
Rn(t)- We observe the same time dependence for both the lengths indicating scaling in the 
late time regime. 

2. Asymmetric case 

In this section, we discuss a ternary mixture with one of the components equally soluble 
in the other two. For this we choose 0\ = cr_i and cr < — (o" 1 + o"_ 1 ) such that the component 
C (i — 0) is the solute. The components A and B are mutually immiscible. We choose A and 
B to have the same volume fraction. We now have A-C and B-C mixtures phase separating 
similar to the binary fluid under critical quench. 

As shown in figure 12, the starting mixed phase separates into percolating domains of 
A-C and B-C phases. Though the early time behavior is similar to the two component 
fluid with a i = <t_i, the late time growth is analogous to the situation where <j\ ^ <r_i 
discussed in section III. The snap shots do not exhibit self-similarity, domains split up and 
combine during coarsening. In this regime, the first zero of the correlation function exhibits 
a dynamic exponent a = 1/2. We would like to point out that for the same density, the two 
component fluid with o"x = <r_i exhibits a late time exponent a = 2/3 and does not show 
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Fig. 13 The first zero of the correlation function (open circles) and the total border length 
(filled squares) in three component asymmetric fluids. Continuous lines with slopes 0.33 and 0.5 

are given as reference 

any evidence of a = 1/2 regime. Thus the third component act to reduce the interfacial 
tension between the A rich and B rich phases leading to fluctuation induced breakup and 
reorganization of the domains. As the volume fraction of the component C is reduced, we 
observe a gradual change from the a=l/2toa = 2/3 regime. 

In order to investigate the validity of the scaling hypothesis in this regime, we compare 
the length scales C(t) and R(t). This is depicted in figure 13. It is evident from the figure 
that these two length scales vary differently at late times. 



V. CONCLUSION 



We have presented the dynamics of domain growth in binary and ternary immiscible 
fluids , using a two-dimensional hydrodynamic lattice gas model. Various dynamic regimes 
are investigated by altering parameters like volume fraction, the interaction strength of the 
individual components and a temperature variable. 

We examine the validity of scaling in two component fluid coarsening by comparing the 
dynamic nature of the first zero of the correlation function and the total border length. While 
scaling holds in the droplet coalescence and the inertial hydrodynamic regimes, violation of 
scaling is observed in the fluctuation dominated regime. 

In symmetric ternary mixtures, we show the existence of a coalescence dominated regime 
with a single length scale exhibiting a dynamical exponent a — 1/2. In asymmetric ternary 
mixtures we use the concentration of the solute as a control parameter. We demonstrate 
that at higher solute concentrations, the coarsening is driven by fluctuations and scaling is 
violated. 
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VI. APPENDIX 



A. Stability of the drop 

In this appendix, we look at the stability of a drop with radius R against perturbations in 
its perimeter. We consider a drop of radius R with the fluid in the region r < R characterized 
by density p, and kinematic viscosity v = - p surrounded by another fluid with density p 
and kinematic viscosity v . For an incompressible and vorticity free fluid we can define the 
velocity component using a streaming potential tjj ss v r = v<j, = — |^ 

Navier-Stokes equation then reduces to 



(d t -uV 2 )V 2 ij(r,<p,t) =0 (17) 
We look for solutions of the former, 0, t) = ipi(r, 0, t) + ijj^r, 0, t) such that 

V 2 ^i = 0. (18) 

(d t -uV 2 )V 2 Mr,<P,t) = 0. (19) 
The general solutions of these equations in cylindrical polar coordinates is 

^(r,0,t) = (a m r m + ^)e^'. (20) 

Mr, 0, t) = (a m I m (kr) + f3 m K m (kr))e m ^e^ . (21) 
Where k = (cu/u) 1 ^ 2 and I m (kr) and K m (kr) are spherical Bessel functions. 



The coefficients a m , b m , a m and (3 m are obtained from the four equations; the continuity 
of v r at r = R , continuity of v<j, at r = R, continuity of tangential stress o r ^ at r = R and 
the relation between normal stresses 



0- rr \ r=R - +p = (Trr\r=R+ + P + + Pa- (22) 

For the pressure p we write P = p( r ) e * m </ > e w *. The pressure due to surface tension is 
obtained from the Young Laplace formula 

„ aim 3 , . 

(r = -^A<R- (23) 

KqOJ 
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These equations can be written in a matrix form, 



R m ~l/R m Ira ~K m 

FT- 1 l/R m+1 -IJR + ilm-i -KJR + £X m _l 

-2(m 2 + m)it— 2 2 -^±^ A 33 A u 

An A12 A13 A14 



(24) 



where, 

A33 
A34 
An 

A12 
A43 
A 44 



2m m 



mk 



~ ( T + ^ + ^ )/m + 2/2 W + ( 2^ + _ kH 

,k 2 2m m 2 N „ mk „ ,m&; k S r . 



k 



m-2 



z/(2m 2 + A; 2 it> 2 - 1 + ^)R m - 2 
Ruirj 

v{2m 2 + k 2 R 2 - 1) 

R(.m+2) 

2m 2 4m 1 am 3 s 
V "~W ~ ~W + R 2 + Rh^' 
, 2m 2 4m 1 T ^ 



(25) 



k 2k 

I m „(1 - 2m)I m -i + -^(m - l)I m -i) 

mrC it 

k 2k 
— (1 - 2m)K m „i + — (m - l)X m+ i) 



In the limit of large viscosity, we can neglect all terms of order \J\ijj/r]). The condition 
for non trivial solutions in this limit leads to a relation between 00 and m; 

m 3 <7 

U = (l-4m 2 )R V (26) 

Since the right hand side of this equation is always negative, we come to the conclusion 
that in the limit of high viscosity there are no unstable modes. In the more general case, 
the condition for non trivial solutions leads to a transcendental equation for u and m. A 
numerical analysis on this equation does not yield any unstable solutions. 
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